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Abstract 

We derive a generalized version of the work fluctuation theorem for nonequilibrium 
systems with spatio-temporal temperature fluctuations. For x 2 -distributed inverse 
temperature we obtain a generalized fluctuation theorem based on g-exponentials, 
whereas for other temperature distributions more complicated formulae arise. Since 
(/-exponentials have a power law decay, the decay rate in this generalized fluctuation 
theorem is much slower than the conventional exponential decay. This implies that 
work fluctuations can be of relevance for the design of micro or nano structures, 
since the work done on the system is relatively much larger than in the conventional 
fluctuation theorem. 



1 Introduction 



In this paper two recent developments are being combined: Fluctuation the- 
orems [1] [13] and Tsallis statistics (and further generalized statistics) in 
non-equilibrium stationary states [14] — [26]. 

The conventional fluctuation theorems give the ratio of the probability P(W T ) 
to find a fluctuation of, e.g., the work W T done on a system over a time interval 
r to the probability P(—W T ) of an equal amount of work —W T done by the 
system. This ratio is given by an exponential function of W T , i.e. one has 
P(W T )/P(— W T ) = exp((3W T ), where (3 is the inverse temperature. On the 
other hand, the 'superstatistical' generalization, derived here, can yield for 
this ratio (for example) a power law instead of exponential behavior. 
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Superstatistics is a 'statistics of a statistics' relevant for driven nonequilibrium 
systems with a stationary state [23] - - [26]. One assumes that the system has 
spatio-temporal temperature fluctuations on a relatively large scale, due to 
external forces acting on the system. One statistics is then given by ordinary 
Boltzmann factors (in space-time cells where the temperature is sufficiently 
constant), and the other one by the probability distribution of (3 in the various 
cells. It has been shown that these types of models generate Tsallis statistics if 
j3 is x 2 -distributed [16,17]. But many other generalized statistics are possible 
as well [22,23,27]. 

We will present a superstatistical generalization of the fluctuation theorem. 
Our result could have important consequences for the work fluctuations in 
micro or nano systems, because the ratio P(—W T )/P(W T ) of the above men- 
tioned probabilities is much larger for the superstatistical case. In fact, this 
implies that for large fluctuations a micro or nano structure will be subjected 
relatively to much more work than in the conventional case. This could, in par- 
ticular, be of relevance if micro or nano systems are placed into a fluctuating 
or turbulent environment, with externally produced temperature fluctuations 
due to constantly acting external driving forces. Moreover, besides the work 
fluctuation theorem, there is also a heat fluctuation theorem [12]. For small 
systems the heat fluctuations allowed by the heat fluctuation theorem can 
also lead to internal fluctuations of the temperature. All this illustrates the 
need to consider fluctuation theorems for systems where the temperature is 
not constant but fluctuates. 

The Fluctuation Theorem can be visualized by considering a very long trajec- 
tory in the phase space of the system in a non-equilibrium stationary state, 
then partitioning this trajectory in segments, on all of which the system spends 
an equal time r. A histogram is then made of W T and W_ T for all segments, 
which will allow an (approximate) determination of the probabilities P(W T ) 
and P(— W T ) and their ratio. 

The main difference between the conventional and the superstatistical Fluc- 
tuation Theorem is that e.g. in the case of a Brownian particle, as considered 
here, the friction and the noise strength in the Langevin equation are con- 
stants, while in the case of superstatistics these quantities can fluctuate in 
space and time. A concrete physical interpretation of the Langevin equation 
with constant coefficients in terms of electric circuits has been proposed in 
[13]. In this approach, the velocity of the Brownian particle corresponds to 
the time-dependent current measured in the circuit. Applying this physical 
interpretation to our Langevin equation with fluctuating coefficients, our su- 
perstatistical ffuctation theorem describes the work fluctuations in a small 
circuit system whose surrounding temperature fluctuates on a relatively large 
time scale. 
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In this paper we will only be concerned with fluctuations of the work done 
on or by a system in a non-equilibrium stationary state in the large time t 
limit. The fluctuations of the heat are much more complicated and will not be 
considered here [12]. There are very many other similar fluctuation theorems 
for work and heat e.g. for transient states and for finite times, all of which 
could be generalized to superstatistical fluctuation theorems. As said, we will 
restrict ourselves here to the simple case, in order to elucidate the principle of 
the superstatistical generalization of the conventional fluctuation theorems on 
the above described model: the replacement of an exponential by e.g. a power 
law. 

This paper is organized as follows. In section 2 we formulate a superstatistical 
Langevin model for the joint process of velocity and position of a Brownian 
particle moving in a changing environment. In section 3 we consider the over- 
damped case of that model, which leads to a superstatistical Langevin equa- 
tion for the position only. In section 4 we derive the superstatistical version 
of the work fluctuation theorem. Finally, in section 5 some concrete examples 
of superstatistical fluctuation theorems are written down. 



2 Superstatistical Langevin model 

Consider a spherical Brownian particle in three dimensions with radius R and 
mass m. We assume that this moves in a fluid of viscosity rj in an environment 
that has spatio-temporal temperature fluctuations. These temperature fluctu- 
ations take place on a rather large spatio-temporal scale and are produced by 
external driving forces that act on the system and change the environment 
of the Brownian particle. A typical example could be a particle moving in 
a turbulent flow experiment. An effective statistical mechanics description of 
the test particle is given by a superstatistics [23], a 'statistics of a statistics', 
where one is given by a probability density associated with the Brownian par- 
ticle for a locally constant temperature, while the other one is given by the 
distribution f((3) of the inverse temperature (3 in the various spatio-temporal 
regions (cells) of approximately constant f3. The temperature [3~ l surround- 
ing the Brownian particle on its way through these spatial regions of size 
varies on the rather large time scale tl (due to both, the spatial movement 
of the particle and also due to a possible explicit time dependendence of the 
temperature field), tl is the typical time scale the Brownian particle spends 
in the cell of size L 3 . We assume that the relaxation time r r of the Brownian 
particle satisfies r r « tl so that local equilibrium can be established. 

To establish a fluctuation theorem, we assume that the motion of the Brownian 
particle is restricted by a time-dependent harmonic potential [9] . For t < 
the minimum x* of the harmonic potential is at the origin, x* — 0, whereas 
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for t > it moves with a velocity vjf. This velocity vj* can, in principle, be an 
arbitrary function of time [11]. The equations of motion for the particle are 



±t = v t (la) 
mvt = -avt-k(xt-x5) + Ct, ( lb ) 

where x t and v t are the position and velocity of the Brownian particle, re- 
spectively. 

There are three different forces acting on the particle. Firstly, the damping 
force -«v t , with Stokes' law yielding a as 

a = 6nr}R. (2) 

Secondly, there is the harmonic force — /c(x t — x£). Finally, there is the random 
force ( t , which is taken to be Gaussian white noise: 

(C*> = 0; (C*&> = 2/rW(* - s) (3) 

The important point is that (3 is not assumed to be constant, but varies on 
the rather large time scale tl. 

The local equilibrium distribution, for some given local /?, assuming that this 
P persists locally for a sufficiently long time interval for the Brownian particle 
to reach local equilibrium, is given by 

WX,V) = (M™) ^(Hvl^lx-x^). (4) 

In the long-term, the particle moves through spatial regions which have differ- 
ent inverse temperatures (3. Denoting the probability density of (3 by f((3), one 
asymptotically obtains a convolution of the various equilibrium distributions 
for different f3. The observed long-term distribution (the marginal distribution 
where (5 is not fixed anymore), relevant for t » t l , is obtained by integrating 
over all (3: 

roc 

p(x,v)= / /(/3)/,^(x,v)d/3. (5) 
Jo 

Let us consider an important example, namely a x 2 -distributed (3 [28,16,17,23]. 
This means the probability density of f3 is given by 

n 

n denotes the number of independent Gaussians X { contributing to the x 2 ~ 
distribution of (3 = Z)" =1 Xf. The mean of the inverse temperature is given by 
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(P) = / Pf(PW = & (7) 
Jo 

and the variance by 

(P 2 ) ~ Pi = ~Pl (8) 

Th 

Let us define the energy of the Brownian particle as 

E= l -m\v\ 2 + l -k\*-*l\ 2 . (9) 

The integration in Eq. (5) can be performed explicitly, and one obtains [17] 
the generalized canonical distributions of nonextensive statistical mechanics 
[14,15,19,20] 

p(x, v) ~ (10) 

(l + ^-l)^)'- 1 

with the following identifications 
and 

P = ^ r. (12) 

l-3(g-l) 1 J 

The marginal distribution (10) describing the long-term distribution of posi- 
tion and momentum of the particle becomes a so-called g-exponential [15,20]. 
For other distributions f(P), one obtains other superstatistics, which all reduce 
to Tsallis statistics for sharply peaked temperature distributions, as shown in 
[23]. 



3 Simplification for the overdamped case 

To proceed to a simpler stochastic model, we may proceed in the usual way 
[29] and consider the system in the overdamped case 

mk < a 2 . (13) 

In that limit one obtains a Langevin equation for the position only, 

i t = -T- 1 (x t -^) + a- 1 C t , (14) 

with relaxation time 



r r = l (15) 
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The local equilibrium distribution reduces to 

Pe^x) = /dv/ e9 (x,v) = (kp/2nf 2 e-^-^\ (16) 
and again the marginal distribution is given by 

p(x) = / f{0) Peq ,p{x)dp. (17) 
For the example of a x 2 -distributed inverse temperature, we obtain 



p(x) ~ -j- (18) 

(i + i^-^jx-x;! 2 )- 1 



with 



and 



^T^I)' (20) 



4 Generalized work fluctuation theorem 

In this section we follow reference [11] with appropriate modifications to solve 
our superstatistical model. 

The work W T done on the system during a time interval r is given by 

Wr = [ T dtv* t -[-k(x t -x* t )}. (21) 
Jo 

The conventional fluctuation theorem for stationary states with a given con- 
stant inverse temperature (3 reads for sufficiently large r: 

p{w t ) • 1 } 

In [11] the fluctuation theorem was derived from a Langevin model with con- 
stant temperature. We now extend this approach to nonequilibrium systems 
with spatio-temporally fluctuating temperatures. 

Let us first restrict to spatial regions (cells) where the temperature can be 
taken to be locally constant. Within these cells, all relevant random vari- 
ables are Gaussians provided r < t l . In Eq. (21), W T is a linear function of 
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x t . Combined with the Gaussian nature both of the Green's function of the 
Ornstein-Uhlenbeck process and of the initial distribution, this implies that 
the distribution P of W T in a given cell for a given time interval r is Gaussian 
as well: 

Here M T denotes the mean of W T and V T the variance of W T . 

Following section C in [11], we obtain straightforwardly for sufficiently large 

T 

M T = k f T dt' 2 f t2 dt[ e-^-'iV^v*, • vi (24) 

JO JO 2 1 

and 

V T = 2(3~ 1 M T . (25) 

This relation proves a local fluctuation theorem in the cells with constant (5. 
Combining Eq. (23) and Eq. (25), we may write 



p ^-^M-iu {w - M ^\ (26) 



and 



p ^ = iik exp {-m { ^ +M ^\' (27) 



which immediately leads to 



P(-W T ) 

P(W T ) 



(2* 



Let us now take into account the temperature fluctuations. We have the fol- 
lowing ordering of time scales 

r r «r < r L « t, (29) 

where r r is the relaxation time to local equilibrium, r is the time scale of a 
local fluctuation theorem, tl is the time scale the particle spends in a cell of 
size L 3 , and t is the total time elapsed. Two different experimental setups seem 
possible, for which we derive two different versions of a generalized fluctuation 
theorem for t — > oo. 



Setup 1. In the experiment one measures the ratio P{— W T ) / P(W T ) in local 
time intervals r < r L , within which the temperature stays sufficiently constant. 
In this case one obtains in the long-term (t » r) 
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Here the expectation (■ • • ) denotes an expectation with respect to the distri- 
bution f(P). 



Setup 2. In a different experimental setup P(W T ) and P{— W T ) are each 
measured over a very long time period that includes many different states with 
different temperatures. In this case one obtains in the long-term (t » r) 

(P(-W T )) f °° f{0) {-^kiWr + M r ) 2 } ygg 

(P(W T )) exp {-^;(W T - M r ) 2 } VW' 



5 Some examples 

5.1 x 2 - distribution 

The assumption of a x 2 -distributed inverse temperature /3 of the form (6) leads 
to a Tsallis-generalized fluctuation theorem. One obtains for an experimental 
setup of type 1 a generalized fluctuation theorem of the form 

, P(-W T ), foo 

{ P(W T 

with 



) = / e-^f(P)dp = (1 + (3(q - l)W T y—, (32) 
J o 



g = 1 + - (33) 

and 

/3 = A). (34) 

This means the ordinary exponential e x in the fluctuation theorem is replaced 
by a g-expontial e x q := (1 — (q — l)x)~~ [19,20]. Note that in the fluctuation 
theorem there is now an asymptotic power law in W T for large W T , which 
means that events with negative work-production, i.e. with work done on the 
system, are now much more likely than in the case without f3 fluctuations. 
Thus rare events are significantly enhanced by the temperature fluctuations: 
previously they were exponentially suppressed, now they obey a power law. 
This can be of potential interest for small electrical circuits [13] in a fluctuating 
environment [23] . Also note that the q relevant for the generalized fluctuation 
theorem as given by eq. (33) is different from the one for the nonextensive 
equilibrium distribution given in Eq. (11) or Eq. (19), due to the fact that 
for the superstatistical invariant distributions we integrate over /5-dependent 
normalization factors, whereas for the fluctuation theorem (in setup 1) we do 
not. 
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For small W T , the result (32) can also be written as the following power ex- 
pansion (cf.[23], eq. (14)): 

<^y> = e-* Wr (i + \(a - i)0X 2 i) 2 /W + •••) (35) 



For setup 2 we obtain the formula 



where 



and 



4M T 



(P(W T )) \(l + p(q - 1)- 



4M T 



(36) 



«= 1+ ^fl (37) 

* = r^fe' (38) 



5.£ Log-normal distribution 



The log-normal distribution 

yields yet another possible superstatistics. m and s are parameters. The av- 
erage (3q of the above log-normal distribution is given by /3q = my/w and the 
variance by a 2 = m 2 w(w — 1), where w := e s [23]. Lognormal superstatis- 
tics has applications in Lagrangian turbulence models [25,26]. Experimentally 
measured data of single-particle accelerations are described very well by these 
models. The physical meaning of the variables occurring in the turbulence 
models is different, for example the variable x in the Langevin equation is 
representing the acceleration of a test particle in the turbulent flow, rather 
than the position of a Brownian particle, while /3 is not related to the inverse 
physical temperature in the flow but to the fluctuating energy dissipation in 
a cascade. Nevertheless, the mathematics is the same and one obtains in ex- 
perimental setup 1 a fluctuation theorem of the form (30), which for small W T 
has the power law expansion (cf. [23], eq.(16)) 

<^y> = + - l)W 2 + 0{W% (40) 

One might speculate that these or similar types of fluctuation theorems de- 
scribe the backward scattering in the turbulent energy cascade. 
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5.3 F- distribution 



Consider a (3 e [0, oo] distributed according to the F-distribution [28,22,23] 
m = n(v + v)/2) fbv\ W2 ^ (41) 

Here w and v are positive integers and b > is a parameter. The average of 
(3 is given by 

fa = w/b(w - 2) (42) 

and the variance by 

a 2 = 2w 2 (v + w- 2)/b 2 v(w - 2) 2 (w - 4). (43) 

For small W T we obtain in experimental setup 1 an expansion of the form 
(cf.[23], eq. (20)) 

<^^> = + + 0{Wl)\. (44) 



5.4 Sharply peaked distributions 



While in general the large W T behavior in the various fluctuation theorems 
depends strongly on the superstatistics considered (i.e. on the function f(/3)), 
the small-WV behavior is universal, in the sense that there is always the same 
quadratic correction term to the ordinary fluctuation theorem. To see this we 
can adopt a proof previously described in [23]. For any distribution f{(3) with 
average (3 := {(3) and variance a 2 := (f3 2 ) — f3 2 we can write 



= e -PoW T ^ e -((3-p )W T ^ 

= e-* w * (l + \a 2 W 2 + f^tSL^p _ p o y )W ^ . (45) 

The coefficients of the powers W[ are the r-th moments of the distribution 
f((3) about the mean (3 . For small enough aW T , in setup 1 the first order 
correction term to the ordinary fluctuation term is always quadratic in W T : 

(^^)^e-^(l + 1 r 2 W 2 ) (46) 
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For sharply peaked distributions in (3, the fluctuation theorem for experimental 
setup 2 also simplifies and one obtains 



(47) 

where (3q and a 2 are the mean and variance of the probability density f{(5) : = 
Cf3 l / 2 f((3) used for mapping type-B superstatistics into type-A superstatistics 
[23]. 
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